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Abstract 



A model-independent search is presented for the production of heavy resonances de- 
caying into top-antitop quark pairs, seen as excesses above the standard model pre- 
diction in the tt invariant mass spectrum. Events containing one lepton (muon or 
electron) and at least two jets are selected from data samples corresponding to 4.4- 
5.0 fb -1 of integrated luminosity collected by the CMS experiment in pp collisions at 
y/s = 7 TeV at the LHC Results are presented from the combination of two dedicated 
searches optimized for boosted production and production at threshold. No excess of 
events is observed over the expected yield from the standard model processes. Top- 
color Z' bosons with narrow (wide) width are excluded at 95% confidence level for 
masses below 1.49 (2.04) TeV and an upper limit of 0.3 (1.3) pb or lower is set on the 
production cross section times branching fraction for resonance masses above 1 TeV. 
Kaluza-Klein excitations of a gluon with masses below 1.82 TeV (at 95% confidence 
level) in the Randall-Sundrum model are also excluded, and an upper limit of 0.7 pb 
or lower is set on the production cross section times branching fraction for resonance 
masses above 1 TeV. 
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1 Introduction 

The top quark is the heaviest known fermion, making it a powerful benchmark to extend our 
understanding of the origin of mass. Because of its large mass, the top quark plays a central 
role in several theories beyond the standard model (SM). These theories predict the existence 
of heavy resonances that manifest themselves as an additional resonant component to the SM 
tt production. Examples of such resonances, which decay preferentially into tt, include mod- 
els with massive color-singlet Z-like bosons in extended gauge theories JlHSl, colorons 0HZI 
or axigluons [8, 9j, models in which a pseudoscalar Higgs boson may couple strongly to top 
quarks ||T0|, and models with extra dimensions, such as Kaluza-Klein (KK) excitations of glu- 
ons |[TT| or gravitons [12J in various extensions of the Randall-Sundrum model ||T3| . 

Recent models HH-HU aimed at explaining the tt charge asymmetry observed at the Teva- 
tron [19-22 J predict resonances in the 0.7-3 TeV mass range with production cross sections of 
the order of a few pb and add renewed interest to the sub-TeV mass region. Independent of 
the exact model, resonant tf production could be visible in the reconstructed invariant mass 
spectrum (Mg). 

Searches performed at the Tevatron have set upper limits on the production cross section of 
narrow resonances (Z' with mass below ~900 GeV) decaying into tt Il23] - [28l . Similarly, searches 
at the Large Hadron Collider (LHC) have set sub-pb limits on the production cross section of 
resonances in the 0.5-3 TeV mass range [|2"9TI3"1"T| . 

In this paper, we present a model-independent search for the production of heavy resonances 
decaying into tt using data collected by the Compact Muon Solenoid (CMS) experiment in 
pp collisions at \/s = 7 TeV at the LHC. Using samples corresponding to an integrated lu- 
minosity of 4.4-5.0 fb -1 , we focus on the semileptonic tt decay mode tt — > (W + b)(W~b) — > 
(qiq 2 b) (t~ V(b) (or charge conjugate) wherein one W boson decays to an electron or muon and 
a neutrino, and the other decays hadronically. The range of 0.5-3 TeV in M t | is covered by the 
combination of two dedicated searches: one optimized for resonances with masses smaller than 

1 TeV (threshold region), and a second one optimized for masses larger than 1 TeV (boosted re- 
gion). Both regions increase the sensitivity of the search by identifying jets originating from the 
hadronization of b quarks (b jets), and separating the samples into various categories depend- 
ing on the lepton flavor, the number of jets, and the number of b jets. The resulting samples are 
dominated by SM tt and W bosons produced in association with jets. A limit on the produc- 
tion cross section of heavy resonances is extracted by performing a template-based statistical 
evaluation of the reconstructed distribution. 

The CMS detector is briefly described in Section |2j Section |3| provides details on the data and 
simulated samples used in the analyses. Sections |4] and [5] describe the event selection and 
the tt event reconstruction, respectively. The main sources of systematic uncertainty in the 
analyses are described in Section |6j Results are shown in Section [7] and a summary is provided 
in Section [U 

2 The CMS Detector 

The central feature of the CMS apparatus is a superconducting solenoid, 13 m in length and 
6 m in diameter, which provides an axial magnetic field of 3.8 T. The bore of the solenoid is 
outfitted with various particle detection systems. Charged particle trajectories are measured 
by the silicon pixel and strip trackers, covering < (p < 2n in azimuth and \rj\ < 2.5, where 
pseudorapidity rj is defined as rj = — ln[tan (8/2)], with 9 being the polar angle of the trajectory 



2 



3 Data and Simulated Samples 



of the particle with respect to the counterclockwise beam direction. A crystal electromagnetic 
calorimeter and a brass/ scintillator hadronic calorimeter surround the tracking volume. In 
this analysis the calorimetry provides high-resolution energy and direction measurements of 
electrons and hadronic jets. Muons are measured in gas-ionization detectors embedded in the 
steel return yoke outside the solenoid. The detector is nearly hermetic, allowing for momentum 
balance measurements in the plane transverse to the beam directions, which are used to infer 
the presence of neutrinos in events. A two-tier trigger system selects the most interesting pp 
collision events for use in physics analysis. A more detailed description of the CMS detector 
can be found in Ref. |32| . 

3 Data and Simulated Samples 

The data analyzed for the threshold analyses were recorded with triggers requiring a single 
isolated (defined in Section |4T] | muon or electron with a transverse momentum (pj) threshold 
of 17GeV or 25GeV, respectively, in combination with a number of jets with a pj threshold of 
30 GeV. Events containing an electron were required to have three or more jets throughout the 
data-taking period, while the minimum number of jets in events containing a muon increased 
from zero to three as the instantaneous luminosity increased. The data analyzed for the boosted 
analyses were recorded with triggers requiring one muon with a pj threshold of 40 GeV or 
one electron with a pj threshold of 65 GeV, with no isolation requirements on the leptons. To 
avoid too high a trigger rate, the electron trigger was prescaled for the highest instantaneous 
luminosities. This resulted in a loss of 0.6 fb _1 of integrated luminosity for the boosted electron 
analysis compared to the other channels. No additional requirements were made on the jets or 
missing transverse energy in the triggers used for the boosted analyses. 

Offline, we use a particle-flow 1133 1 based event reconstruction, which combines information 
from each subdetector, including charged particle tracks from the tracking system and de- 
posited energy from the electromagnetic and hadronic calorimeters, to reconstruct all parti- 
cles in the event. Particles are classified as electrons, muons, photons, charged hadrons, and 
neutral hadrons. Particles identified as originating from multiple primary collisions at high 
instantaneous luminosity (pileup) are removed from the event. 

Muons are reconstructed using the information from the muon chambers and the tracking de- 
tectors |34|. Tracks are required to have at least 11 hits including at least one in the pixel layers. 
The tracks must also pass within 0.02 cm of the beam spot in the plane transverse to the beam, 
and within 1 cm along the beam axis. 

Electron candidates are initially identified by matching a track to a cluster of energy in the 
electromagnetic calorimeter. Candidates are selected [35] using shower-shape information, the 
quality of the track and the spatial match between the track and electromagnetic cluster, the 
fraction of total cluster energy in the hadronic calorimeter, and the amount of activity in the 
surrounding regions of the tracker and calorimeters. Electrons coming from photon conver- 
sions in the detector material are rejected if there are missing hits in the inner tracker layers or 
if there is another close track with opposite charge and with a similar polar angle. 

Jets are reconstructed by clustering the particle-flow candidates not identified as leptons using 
an anti-fcx algorithm with a distance parameter R = 0.5 11361 . Corrections are applied to account 
for the dependence of the detector response to jets as a function of tj and pj fl37| and the effects 
of pileup. The jets associated to b quarks are identified using an algorithm that reconstructs 
the secondary vertex corresponding to the decay of a B hadron. When no secondary vertex 
is found, the significance of the impact parameter with respect to the primary vertex of the 
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second most displaced track is used as a discriminator to distinguish decay products of a B 
hadron from prompt tracks l|38| . 

The negative of the vector sum of the momenta of all reconstructed particles in the plane trans- 
verse to the beam is the missing transverse momentum f39l , with magnitude denoted by miss- 
ing transverse energy E™ 1SS . 

The SM background processes are simulated by MadGraph 5.1.1 001, pythia 6.4.24 EH, and 
POWHEG [42] event generators using CTEQ6L parton distribution functions of the proton [43J. 
The generated events are subsequently processed with PYTHIA to provide the showering of the 
partons and fully simulated with CMS software based on Geant4 |44ll45| . 

The W boson and Drell-Yan production in association with up to four jets are simulated with 
MADGRAPH, with additional jet production described via matrix elements matched to parton 
showers using the MLM prescription |46l with a matching threshold of 20 GeV. The next-to- 
next-to-leading order (NNLO) production cross sections times branching fractions into leptons 
(electrons, muons and taus) are used [47J: 31.3 nb for W, and 3.05 nb for Drell-Yan production of 
dilepton final states with invariant mass > 50 GeV. The background from Drell-Yan production 
of dilepton final states with invariant mass < 50 GeV is negligible. The contribution from QCD 
multijet processes is obtained directly from data. 

The SM tt events are generated with MadGraph, assuming a top-quark mass of 172.5 GeV. 
Higher-order gluon and quark production is described by the matrix elements with up to three 
extra partons beyond the tt system. The chosen threshold for the matching is 40 GeV, which 
ensures a smooth transition from the matrix element to the parton showering description. An 
additional tt sample is generated using POWHEG to provide a cross-check and to estimate 
systematic uncertainties in the modeling. The inclusive tt cross section value of 157.5 pb is 
used BUSS). 

Single top-quark production is modeled in POWHEG. The approximate NNLO cross sections of 
42 pb and 3.2 pb are used for i-channel and s-channel single top-quark production, respectively, 
along with the corresponding single t-quark production cross sections of 23 pb and 1.4 pb. The 
approximate NNLO value of 7.9 pb is used for Wt and Wt associated production I5Q]452|. 

Finally, as reference models for new physics, we use the sequential standard model (SSM) top- 
color Z' bosons with a natural width Yjj equal to 1.2% (narrow width) and 10% of the Z' mass 
mjj based on EHZ1 and KK gluons based on |TT| . Signal samples are generated with PYTHIA 
8.145 with a range of masses between 500 GeV and 3 TeV. Only decays into tt are simulated in 
the Z' samples. The KK gluons are simulated with branching fractions to tf of 0.93, 0.92, 0.90, 
and 0.87 for resonance masses of 1, 1.5, 2 and 3 TeV. 

4 Event Selection 

To study the range of 0.5-3 TeV in M t j, two complementary strategies are pursued: firstly, the 
threshold search focuses on the 0.5-1 TeV mass range using criteria optimized to identify top 
quarks produced with a small boost in the detector frame and hence with well-separated decay 
products. In this region, if all decay products are reconstructed within the kinematic accep- 
tance, we expect the final state to contain exactly one isolated lepton, four jets produced by the 
four quarks (two of which are b jets) in the semileptonic tt decay, and E™ 1SS . 

Secondly, for resonance masses above 1 TeV, the highly Lorentz-boosted top quarks will yield 
collimated decay products that are partially or fully merged. This can be seen in Fig. [TJ which 
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Figure 1: The distribution of the minimum AR of all three possible pairings between the three 
quarks (qi,q2,b) of the hadronic top-quark decay for SM tt production and two different Z' 
mass hypotheses. For events with ARmin smaller than the parameter R = 0.5 in the jet cluster- 
ing, jets merge and fewer than three jets are reconstructed. 

shows that in the boosted region the angular distance between the partons is smaller than the 
jet clustering distance parameter. As a consequence, the products of the hadronically decay- 
ing top quark might be reconstructed as fewer than three jets, and the leptons might not be 
isolated. The boosted search thus selects events containing one electron or muon with no iso- 
lation requirement and at least two jets. 

4.1 Threshold analyses 

We select events containing either one isolated muon with pj > 20GeV and \tj\ < 2.1, or one 
isolated electron with pj > 30GeV and \t]\ < 2.5. The isolation requirement is based on the 
ratio of the total transverse energy observed from all hadrons and photons in a cone of size 

AR = y (A(p) 2 + (At]) 2 < 0.4 around the lepton direction to the transverse momentum of the 
lepton itself. This quantity is required to be less than 0.125 for muons and less than 0.1 for 
electrons. Events with two isolated lepton candidates are vetoed to reduce the background 
from Drell-Yan and tt production in which both W bosons decay leptonically. 

Events are further required to contain at least three jets with \t]\ < 2.4 and pj > 50GeV, and 
additional jets with \rj\ < 2.4 and pj > 30GeV, if any. To enhance the rejection of background 
from W-boson and Drell-Yan production in association with relatively low-px jets, the leading 
jet is required to have pj > 70 GeV. Multijet background is suppressed further by requiring 
£miss > 20 GeV. The fraction of simulated semileptonic signal events passing this selection 
varies from 16 to 35% for resonance masses below 1 TeV. 

Events are then separated into eight categories according to the lepton flavor (electron or 
muon), the number of jets, and the number of b-tagged jets. The categories defined by jets 
are: events with three jets, of which at least one is b tagged; events with four or more jets, of 
which none is b tagged; events with four or more jets, of which exactly one is b tagged; and 
events with four or more jets, of which at least two are b tagged. 

4.2 Boosted analyses 

We select events containing either one muon with pj > 42 GeV and \t]\ < 2.1, or one electron 
with pj > 70 GeV and \rj\ < 2.5, and at least two jets with \ t]\ < 2.4 and pj > 50 GeV. The lead- 
ing jet pj lower threshold is set to 250 GeV (150 GeV) in the muon (electron) channel. No isola- 
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tion requirement is applied to the leptons. Multijet background is reduced with a requirement 
on the AR separation in the 2D plane: AR (lepton, closest jet) > 0.5 or p^ 1 (lepton, closest jet) > 
25 GeV. Here, p is defined as the magnitude of the lepton momentum orthogonal to the clos- 
est jet axis, where any jet with p T > 25 GeV is considered. We also require the scalar quantity 

L T > 150 GeV, where L T = Ef iss + p l j Pton . 

In the electron channel only, the multijet background is further reduced by requiring that 
£rmss > 50 GeV and applying a series of topological requirements that ensure the missing 
transverse momentum does not point along the transverse direction of the electron (e) or of 
the leading jet (j): 

^— E? iss + 1.5 < A(p{(e or i),Ef ss \ < — E^ ss + 1.5. 

75GeV T ^ u u 1 s 75GeV T 

Even though the lepton pj requirements are dictated by the trigger threshold, the leading jet 
pi requirement is chosen so that the total transverse energy of the event (including E miss ) is as 
close as possible in both channels. In addition, we ensure the two channels contain no overlap 
with each other by vetoing events that contain a second lepton. 

Events are separated into four categories according to the lepton flavor (electron or muon) and 
the number of b-tagged jets: either no b-tagged jets, or at least one b-tagged jet. The fraction of 
simulated semileptonic signal events passing this selection varies from 13 to 24% for resonance 
masses between 1 and 3 TeV. 



5 The tt Event Reconstruction 

The four-vectors of the top quark and antiquark candidates are reconstructed by assigning the 
final state objects in each event to either the leptonic or the hadronic leg of the tt pair decay. We 
then choose between the possible hypotheses using the criteria described below that depend 
on the number of reconstructed jets. This tt reconstruction process results in a unique value for 
the reconstructed for each event. 

First, the charged lepton and the E™ 1SS are assigned to the leptonic leg, where E miss is inter- 
preted as the transverse component of the momentum of the neutrino. Imposing the condition 
that the invariant mass of the lepton and neutrino is equal to the mass of the W boson (80.4 GeV) 
allows the construction of a quadratic equation for the longitudinal component of the momen- 
tum of the neutrino. In the absence of a real solution, the boosted analyses retain the real part 
of the complex solution. The threshold analyses modify the components of E™ ss by the min- 
imal amount in |AE™ 1SS J + |AE^ lss y | to give one real solution, which results in an improved 
mass resolution. If there are two real solutions, hypotheses are built for both cases, effectively 
doubling the number of combinations for that event. 

For events with four or more jets in the threshold analyses, the choices of neutrino solution 
and jet association are made simultaneously by forming a x 2 horn the sum of the normalized 
squared deviations of the leptonic top-quark mass, hadronic top-quark mass, hadronic W mass, 
pi of the tf system, and the ratio of the pi of the four selected jets to the pi of all jets in the event. 
The central values and widths used are obtained from the distributions of these quantities in 
the Monte Carlo simulation. The x 2 is calculated for each possible combination, including the 
two neutrino solutions if they are both physical. The b-tagged jets may only be associated to 
a b quark in the decay chain, thereby reducing the number of possible combinations. For each 
event, the combination with the smallest value of x 2 is chosen. The association of jets to the W 
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5 The tt Event Reconstruction 



boson and the b quarks is found to be correct in approximately 80% of simulated tt events for 
which the four jets in the decay chain are reconstructed. 

For events with only three jets in the threshold analyses, it is assumed that two jets from the 
tf decay may have merged. The leptonic W boson is first reconstructed as described above. 
The solution for the longitudinal neutrino momentum is chosen to give the closest match to the 
leptonic top-quark mass when the leptonic W boson is combined with any of the three jets. The 
invariant mass of the leptonic W and all three jets together is then taken as an estimate of M t j. 

For the boosted analyses, we allow for collimated decay products that are partially or fully 
merged by considering all hypotheses that have exactly one jet assigned to the leptonic leg, 
and at least one jet assigned to the hadronic leg. A two-term x 2 is constructed from the sum of 
the normalized squared deviations of the leptonic top-quark mass and the hadronic top-quark 
mass. For each event, the combination with the smallest value of x 1 (labeled Xmm) i s chosen. 



Next, the event selection described in Section 4.2 is extended by applying additional conditions 



that improve the overall sensitivity of the boosted analyses. For the electron channel only, the 
transverse momentum of the reconstructed leptonic top quark is required to be greater than 
100 GeV. We require Xmm < & for both channels. This value is chosen such that the efficiency 
for this cut is 50% for signal and approximately 10% for the W+jets background. Finally, we 
categorize events according to the number of b-tagged jets as either with no b-tagged jets, or 
with at least one b-tagged jet. 

For the threshold analyses, the multijet background contribution in each of the categories is de- 
termined by defining two control regions in data. The first contains events with the same lepton 
isolation requirement as the signal (less than 0.1) but with the E™ ss requirement inverted (i.e. 
less than 20 GeV). The second multijet-dominated region contains events with either a muon 
with the isolation value between 0.2 and 0.5 and no £™ 1SS requirement, or an electron candi- 
date that has been rejected as consistent with being a photon conversion. The distribution of 
a chosen variable in the multijet-dominated region is used to define a template. This template 
is fitted to the distribution in the low-£^ lss region together with a template for the distribution 
of other backgrounds taken from simulation. The chosen variable is selected based on both its 
power to discriminate between the contributions and the reliability of its modeling in simula- 
tion. For muon events the pj of the vector sum of the jet momenta is used, whereas for electron 
events the reconstructed electron tj is used. This provides the multijet background yield in the 
low-£^ lss control region, which is then used to extrapolate into the signal region with the rel- 
ative yields below and above £ip iss of 20 GeV in the second control region. The second control 
region is also used to define the shape of the multijet background contribution to each cat- 
egory. In the boosted analyses, the multijet contamination after the final selection is found to 
be negligible. 

The numbers of expected and observed events in each analysis channel for the threshold and 
boosted analyses are summarized in Tables [l] and [2j respectively. The Z' samples are normal- 
ized arbitrarily to cross sections times branching fractions of 1 pb. For the threshold analyses, 
the simulated samples are normalized to theoretical predictions. For the boosted analyses, the 
yields of the simulated samples are normalized to data using scale factors derived in a max- 
imum likelihood fit to the M t j distribution in both channels simultaneously, as described in 
Section |7| Figures [2] and [3] show the distributions for the threshold and boosted analyses, 
respectively. Figure [3] also shows the distribution of the number of jets in the events for the 
boosted analyses. It can be observed that, in the boosted region, the signal populates the 2-jet 
bin while the SM background has larger jet multiplicity. Good agreement is observed in all 
cases between data and the SM predictions. 
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Table 1: Number of expected and observed events in the threshold analyses for an integrated 
luminosity of 5.0 fb -1 . The narrow-width 71 samples are normalized to cross sections times 
branching fractions of 1 pb. The other simulated samples are normalized to theoretical predic- 
tions. The uncertainty in the total background corresponds to yield changes originating from 
the systematic uncertainties associated with the jet energy corrections, jet energy resolutions, b 
tagging, and pileup. The statistical uncertainties for the simulated samples are negligible. 



Threshold analyses, muon channel 


Sample 


Nj et = 3 


Njet > 4 


Nj e t > 4 


Njet > 4 




N b .t ag > 1 


Nb-tag = 


N b - tag = 1 


Nb-tag > 2 


71 (M=0.5TeV) 


48.5 


14.0 


41.1 


34.6 


71 (M=1.0TeV) 


68.5 


36.1 


95.5 


74.7 


71 (M=1.5TeV) 


56.4 


33.9 


76.5 


50.9 


71 (M=2.0TeV) 


38.0 


32.4 


60.7 


37.5 


tt 


5612 


2988 


7802 


6093 


W/Z+jets 


1727 


7705 


1296 


173 


Single top 


550 


202 


423 


228 


Multijet 


164 


195 


104 


152 


Total background 


8052 ±511 


11089 ± 1241 


9626 ± 822 


6646 ± 687 


Data 


8465 


10714 


9664 


6697 


Threshold analyses, electron channel 


Sample 


Njet = 3 


Njet > 4 


Njet > 4 


Njet > 4 




Nb-tag > 1 


Nb-tag = 


Nb- tag = 1 


Nb-tag > 2 


71 (M=0.5TeV) 


34.7 


10.5 


29.4 


25.0 


71 (M=1.0TeV) 


58.9 


32.4 


85.0 


67.4 


71 (M=1.5TeV) 


51.2 


31.7 


73.5 


50.8 


71 (M=2.0TeV) 


33.8 


30.2 


59.5 


37.6 


tt 


4307 


2395 


6183 


4770 


W/Z+jets 


1372 


6355 


1051 


142 


Single top 


428 


158 


345 


184 


Multijet 


491 


1398 


504 


210 


Total background 


6597 ± 442 


10307 ±1136 


8083 ± 721 


5306 ± 514 


Data 


6932 


10008 


7946 


5309 



8 



6 Systematic Uncertainties 



Table 2: Number of expected and observed events in the boosted analyses for an integrated 
luminosity of 4.4-5.0 fb" 1 . The narrow-width 71 samples are normalized to cross sections times 
branching fractions of 1 pb. The other simulated samples are normalized to data, see Section 7. 
The uncertainty in the total background corresponds to yield changes originating from the 
systematic uncertainties associated with the jet energy corrections, jet energy resolutions, b 
tagging, and pileup. The statistical uncertainties for the simulated samples are negligible. 



Boosted analyses 


Electron channel 


Muon channel 


Sample 


Nb-tag = 


Nb-tag > 1 


Nb-tag = 


Nb-tag > 1 


71 (M= 1 TeV) 


17.1 


36.5 


27.8 


48.3 


71 (M= 1.5 TeV) 


44.7 


55.4 


95.9 


94.4 


71 (M= 2 TeV) 


62.1 


52.8 


146.3 


94.1 


71 (M= 3 TeV) 


57.2 


36.9 


155.2 


69.0 


tt 


172 


336 


157 


262 


W/Z+jets 


95 


6 


149 


9 


Single top 


9.3 


15 


8.1 


11 


Total background 


276 ± 58 


357 ± 50 


314 ± 72 


282 ± 34 


Data 


277 


354 


300 


269 



6 Systematic Uncertainties 

Systematic uncertainties enter the analyses in two ways: those related to the total normalization 
of the simulated samples, and those from the effects that change both the normalization and 
shape of the background and expected signal distributions. 

Normalization uncertainties on the theoretical production cross sections are considered for all 
background processes. In some instances, larger uncertainties are used for the boosted analy- 
ses as they probe a limited region of phase space. The following variations on the rates, which 
were obtained in a previous analysis |53|, are included: tt (15%); single top-quark for threshold 
(30%) and for boosted (50%) analyses; W/ Z+light-quark jets correlated (50%) and additional 
Drell-Yan uncorrelated (30%) for threshold analyses, W+light-quark jets (50%) and uncorre- 
cted Z+light-quark jets (100%) for boosted analyses; W/Z+heavy-quark jets (100%). In addi- 
tion, a 2.2% uncertainty in the luminosity [|54| and 3% (5%) lepton trigger and identification 
uncertainty is applied to all simulated samples for the threshold (boosted) analyses. 

Several sources of systematic uncertainty affect both the shape and the rate of the templates 
used in the analyses. The uncertainty on the energy of jets is of the order of a few percent and is 
parametrized as a function of the jet pj and rj [37). The uncertainty on the jet energy resolution 
varies from 6 to 20% depending on the jet rj. The effect of both variations is propagated to 
the event E™ 1SS . The uncertainty on the b-tagging efficiency for b jets ranges from 1.6 to 8% 
depending on the jet rj and is doubled for b jets with pj > 670 GeV II38II . The uncertainty on the 
b-tagging efficiency for c jets is taken as twice the uncertainty for b jets. The uncertainty for all 
other jets (mistag rate) is 11%. 

Some of the theoretical uncertainties affect the normalization and shape of the simulated sam- 
ples. A simultaneous variation of the factorization and renormalization scales to half and twice 
the nominal scales is allowed for the tt and W+jets samples. The matrix element to parton 
shower matching threshold and the amount of initial- and final-state radiation are also allowed 
to vary for these samples. A further uncertainty is included as the difference between the tt 
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Figure 2: Comparison of the reconstructed M t j in data and SM predictions for the threshold 
analysis with (a) 3 jets of which > 1 b tagged, (b) 4 jets, none of which is b tagged, (c) 4 jets of 
which one is b tagged, (d) 4 jets of which >2 are b tagged. Expected signal contributions for 
narrow-width topcolor Z' models at different masses are also shown. For clarity, a cross section 
times branching fraction of 20 pb is used for the normalization of the Z' samples. 
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6 Systematic Uncertainties 



CMS, 4.4-5.0 ft)" 1 , Vs = 7 TeV CMS, 4.4-5.0 ft)"', Vs = 7 TeV 




Figure 3: Comparison of the reconstructed M t j in data and SM predictions for the boosted 
analysis with (a) no b-tagged jets, (b) >1 b-tagged jets. Comparison of the jet multiplicity 
distribution in data and SM background predictions for the boosted analysis with (c) no b- 
tagged jets, (d) >1 b-tagged jets. Expected signal contributions for narrow- width topcolor Z' 
models at different masses are also shown. A cross section times branching fraction of 1.0 pb is 
used for the normalization of the Z' samples. 
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production models in POWHEG and MADGRAPH. For all simulated samples, the minimum 
bias cross section is varied by 1 standard deviation of its measured value to account for the 
effect of pileup. 



7 Results 

The statistical analysis is based on a binned likelihood of the distributions in the consid- 
ered channels, i.e., eight channels for the threshold analyses and four channels for the boosted 
analyses. The number of events in bin i is assumed to follow a Poisson distribution with mean 
A„ given by the sum over all considered background processes and the Z' signal. The signal is 
scaled with a signal strength modifier ji, which is the signal cross section in pb: 

\i = }iSi + Y^B kri . 

k 

Here, k runs over all considered background processes, is the background template for back- 
ground k, and S is the signal template, scaled according to luminosity and a signal cross section 
oflpb. 

The presence of systematic uncertainties affects the yields A,-. A nuisance parameter 6 U is thus 
introduced for each independent source of systematic uncertainty considered. A rate-only un- 
certainty is modeled with a coefficient for the template B/ c with a log-normal prior. A rate and 
shape uncertainty is modeled by choosing a Gaussian prior for 8 U and using this parameter 
to interpolate between the nominal template and the shifted templates obtained by applying 
a ±1 a systematic shift to the simulated samples. This interpolation uses a smooth function, 
which is cubic in the range ±1 a and linear beyond ±1 c. 

To set exclusion limits on a Z' hypothesis, we define a test statistic q^, which depends on the 
hypothesized signal rate y. ||55| : 



q ¥ = -2 In 



£(0 /( ,f/|data) 
£(6>,/i|data) ' 



where 8^, fi, and 9 are the values of the parameters 8 and }i that maximize the likelihoods in the 
numerator and denominator, and the subscript y. in 6^ indicates that the maximization of the 
likelihood in the numerator is done under the hypothesis of a signal of strength }i. One-sided 
limits on the Z' production rate are obtained by constraining \i < ]i in the denominator. 

We use the modified frequentist construction CL S |56ll57l to calculate the 95% confidence level 
(CL) upper limits on the Z' — > tt cross section. We define two tail probabilities associated 
with the observed data; the probability to obtain a value for the test statistic q^ larger than the 
observed value ^° bs for the signal+background, and for the background-only hypotheses: 

CL s+h = P( qil >qf s \e rF ) 
CL b = P(q pi >qf s \9, j i = 0) 



and obtain CL S from the ratio 

r , _ CL s+b 

To quote the upper limit on }i at the 95% confidence level, we adjust }i until we reach CL S = 0.05. 
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7 Results 



The expected upper limits are calculated using background-only pseudo-experiments {}i = 0) 
and calculating the upper limit for each pseudo-experiment. The expected limit is given by the 
median of the distribution of upper limits, and the central 68% and 95% give the ±1 and ±2 
standard deviation (s.d.) excursions. 

The number of simulated background events in the M t i > 2 TeV region that pass the boosted 
selection is rather limited. To ensure a proper background modeling in the entire range, we 
merge bins in the M t f distribution requiring a minimum number of background events per bin. 
The bins are chosen such that the uncertainty on the number of expected background events 
due to the limited number of simulated events is not worse than 30% in all channels. The 
uncertainty due to finite size of the simulated samples is taken into account using the "Barlow- 
Beeston lite" method |58] that defines one additional nuisance parameter with a Gaussian dis- 
tribution for each bin, and performs the maximization of the likelihood with respect to these 
new parameters analytically. 

Figures |4] and [5] show the expected and observed 95% CL upper limits for the product of the 
production cross section times branching fraction of hypothesized resonances that decay into 
tf as a function of the invariant mass of the resonance. The dashed lines indicate the values 
predicted by various models for new physics processes. The expected mass exclusion region for 
a topcolor 71 with Y z , fm z , = 1.2% is M z , < 1.53 TeV, the observed exclusion is M z , < 1.49 TeV. 
For wide resonances with T Z i / m Z i = 10%, the exclusion mass region is M z * < 2.04 TeV for both 
the expected and observed limits. In Fig. |4j the vertical dashed line indicates the transition 
between the threshold and the boosted analyses, chosen based on the sensitivity of the expected 
limit. For a Kaluza-Klein excitation of a gluon (gKK) the exclusion mass region is M(gKK) < 
1.82 TeV for both the expected and observed limits. 
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Figure 4: The 95% CL upper limits on the product of the production cross section dj] and the 
branching fraction B of hypothesized resonances that decay into tf as a function of the invariant 
mass of the resonance. The Z' production with r Z //m Z / = 1.2% (a) and 10% (b) compared to 
predictions based on [5|. The ±1 and ±2 s.d. excursions from the expected limits are also 
shown. The vertical dashed line indicates the transition between the threshold and the boosted 
analyses, chosen based on the sensitivity of the expected limit. 
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7 Results 
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Figure 5: The 95% CL upper limits on the product of the production cross section Okk and the 
branching fraction B of Kaluza-Klein excitation of gluon production from fTTH , compared to the 
theoretical prediction of that model. The ±1 and ±2 s.d. excursions from the expected limits 
are also shown. 
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8 Summary 

Results from a model-independent search for the production of heavy resonances decaying 
into tf are presented. The data sample corresponds to an integrated luminosity of 4.4-5.0 fb _1 
recorded in 2011 by the CMS detector in pp collisions at a/s = 7TeV at the LHC. After analyz- 
ing events with a lepton (muon or electron) plus jets final state, no evidence of such massive 
resonances is found above the SM prediction. Therefore, limits are set on the production of 
non-SM particles. Topcolor Z' bosons with a width of 1.2 (10)% of the Z' mass are excluded at 
95% CL for masses below 1.49 (2.04) TeV; an upper limit of 0.3 (1.3) pb is set on the production 
cross section times branching fraction for a resonance mass of 1 TeV. In addition, Kaluza-Klein 
excitations of a gluon with masses below 1.82 TeV (at 95% CL) in the Randall-Sundrum model 
are excluded; an upper limit of 0.7 pb is set on the production cross section times branching 
fraction for a resonance mass of 1 TeV. In both instances, the upper limits are lower for larger 
resonance masses. These results set the most stringent limits, to date, for tt resonant production 
in the 0.5-2 TeV mass range. 
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